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1. Afbakening van het gebied 


Het woord plasma werd in het begin van de dertiger 
jaren ingevoerd door I. Langmuir in zijn beschrijving van de 
positieve kolom of zuil van een glim- of boogontlading. De 
zuil onderscheidt zich van de elektrodegebieden door een vrij- 
wel volledige wederzijdse compensatie van de ruimteladingen 
van de ionen en de elektronen en door de afwezigheid erin van 
sterke elektrische velden. 

Het is gebleken dat theorieën waarmee eigenschappen 
van geïoniseerd gas beschreven worden ook tot op zekere hoogte 
toepasbaar zijn op geleidende vloeistoffen en vaste stoffen. 
Dit leidt ertoe dat men het begrip plasma ook wel eens in een 
ruime betekenis opvat, die al deze toepassingen omvat. In deze 
voordrachten zullen we enkele voorbeelden hiervan aanstippen, 
maar ons overigens beperken tot het plasma in zijn oorspronke- 
lijke betekenis van geïoniseerd gas. | | 

B 
2. Theotetische modellen 


Als men in de mechanica van de deformeerbare media 
de mogelijkheid van een elektrische stroom in het fluïdum in- 
_troduceert, kan er in de krachtenvergelijking een term 5 x B 
optreden, die de wisselwerking met een magnetisch veld be- 
schrijft. Men beweegt zich dan op het terrein van de magneto=- 
hydrodynamica en wel in eerste instantie met behulp van een 


êên-vloeistof theorie, waarin de atomaire structuur van de 
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vloeistof niet aan de orde komt. Deze theorie is evenzeer van 
toepassing op het gedrag van vloeibare geleiders als op dat 
van gasvormige plasma's. 

Een verfijning hiervan houdt rekening met het feit 
dat een geleidend gas in het algemeen bestaat uit drie compo- 
nenten, namelijk elektronen, ionen en neutralen. Wil men dan 
in detail nagaan hoe b.v. een uitwendige kracht als 5 x B op 
het gas werkt, dan kan in een zekere situatie het antwoord 
als volgt luiden. De stroming van de elektronen geeft aanlei- 
ding tot een Lorentz-kracht, waardoor hun snelheid zijdelings 
wordt afgebogen. Dit geeft een ladingsscheiding en dus een in- 
wendig elektrisch veld waardoor de ionen met de elektronen wor- 
den meegesleept. De wrijving tussen ionen en neutralen kan er 
dan tenslotte voor zorgen dat ook de laatste worden meegesleept. 
Dergelijke beschouwingen passen in een drie-vloeistoffen theo- 
rie. Voor volledig geïoniseerde gassen met één soort ionen kan 
men uiteraard volstaan met een twee-vloeistoffen theorie. 

Vooral bij zeer langzame oscillaties in plasma's met 
zeer hoge botsingsfrequenties is de koppeling tussen de compo- 
nenten zo sterk dat ze elkaar in hoge mate kunnen meeslepen; 
dan is het gebruik van de EEn-vloeistof magnetohydrodynamica 
verantwoord. Aan de andere kant zullen bij zeer snelle oscilla- 
ties de zware deeltjes door hun grote traagheid nauwelijks in 
beweging komen. Dergelijke verschijnselen kan men goed beschrij- 
ven door de ionen als stilstaand te beschouwen; deze fungeren 
dan alleen om de ruimtelading van de elektronen in de even- 
wichtstoestand te neutraliseren. Dit model is ook van toepas- 
sing op metalen en halfgeleiders. Men noemt daarom ook de ge- 
leidingselektronen in vaste stoffen wel eens een plasma. 

Als het gaat om een beschrijving van in het labora- 
torium of in het heelal optredende verschijnselen zal men in 
het algemeen geen theorie kiezen die meer details aankan dan 
die waarin men geïnteresseerd is, daar de mathematische compli 
caties al spoedig bezwaarlijk worden. Aan de hierboven genoem- 
de continuüm-theorieën heeft men ook veelal wel zijn handen 
vol. Er bestaan echter theoretische finesses waarvoor ze on- 


toereikend zijn, en waarvoor men de kinetische gastheorie te 
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hulp moet roepen. Deze beschrijft het plasma in termen van de 
snelheidsverdeling als functie van plaats en tijd. Door de des- 
betreffende grootheden te vermenigvuldigen met de verdelings=- 
functie en over de snelheid te integreren, krijgt men b.v. de 
ladingsdichtheid en de elektrische stroomdichtheid. De resul= 
terende inwendige velden verschijnen samen met de uitwendige 
velden in de Boltzmann-vergelijking die zegt hoe de verdelings- 
functie onder invloed van deze velden verandert. We zullen ons 
echter in deze voordrachten beperken tot enkele voorbeelden 


van de continuüm-beschrijving van het plasma. 


3. Magnetohydrodynamica (één-vloeistof model) 
3.1 Gevangen Flux 


Het plasma wordt als elektrisch neutraal beschouwd, 
zodat er in de bewegingsvergelijking geen elektrische kracht 
optreedt. Wel kunnen er stromen lopen, waardoor er magnetische 
krachten van de gedaante 5 x B kunnen optreden. Veelal wordt 
viscositeit verwaarloosd en wordt de druk isotroop verondersteld. 
In dit hoofdstuk komt de vraag of het plasma een compressibele 
of incompressibele vloeistof is niet expliciet aan de orde; de 
dichtheid wordt als constant verondersteld, en er worden geen pro- 
cessen beschreven die in een compressibel medium aanleiding zouden 
geven tot div v # O. Het plasma is volledig gekarakteriseerd door 


op en n. De Navier-Stokes-vergelijking gaat dan over in 


=j X B -— Vp (1) 
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(In astrofysische toepassingen heeft men rechts ook nog een term 
e ovy ten gevolge vàn de zwaartekracht). 
De stroomdichtheid hangt samen met het elektrische veld 
volgens | 
E + vxBe= nj (2) 


waarin het linkerlid beschouwd kan worden als het elektrische veld 
dat wordt waargenomen in het coördinatenstelsel dat met snelheid v 
beweegt, terwijl E en B in een stilstaand stelsel gemeten zijn. 
Het stelsel vergelijkingen wordt gecompleteerd door de Maxwell- 
vergelijkingen, de continuïiteitsvergelijking en de toestandverge- 
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lijking voor het gas of de vloeistof. 


> f 


EET NE EP MOE OT De ENE FE a: ME ad: od ew SER ad ne Oben OE OOB en 


In de Maxwell-vergelijkingen wordt in het algemeen u = u ene =e 


©) 


gesteld, wat er op neerkomt dat kringstromen en ladingsscheiding 


niet verdisconteerd worden in u en e,‚ maar expliciet verschijnen 


als Ef (We zullen echter een voorbeeld ontmoeten waarin dan toch 
achteraf de termen j en 9É/3t kunnen worden samengetrokken door een 


EF Én te introduceren.) 


In de "ideale" 


n = 0, waardoor volgens 


magnetohydrodynamica stelt men ook nog 


(2): 


B=0 … (3) 


Beschouw nu een gesloten contour met oppervlakte O die overal 


meebeweegt met de lokale vloeistofsnelheid. De magnetische 


flux door deze contour zij ®. Dan geldt (zie figuur 1): 
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“Gevangen flux”. 


waarin de eerste term rechts veroorzaakt wordt door de veran= 
dering van de magnetische veldsterkte en de tweede door de be- 
weging van de contour. De kringintegraal wordt volgens Stokes 


getransformeerd tot een oppervlakte-integraal, terwijl volgens 


Maxwell SRS == Vv X E, zodat met (3): 
mm 
En zelf Wx È + vx B|.aó = 0. (4) 
Opp. 


De magnetische flux door onze contour is derhalve constant; men 
spreekt van gevangen of —- merkwaardigerwijs — van ingevroren flux, 
en, omgekeerd, van buitengesloten flux. Door een magnetische kracht- 
lijn kan men twee elkaar in die krachtlijn snijdende magnetische. 
oppervlakken leggen, dat zijn vlakken die worden gevormd door kracht- 
lijnen. Een willekeurige contour in één van die vlakken omsluit geen 
flux. Bewegen die vlakken met de vloeistof mee dan blijven ze flux- 
vrij, m.a.w. ze blijven magnetische oppervlakken en hun snijlijn | 
blijft een krachtlijn. Alle vloeistofelementen op één krachtlijn 
blijven dan steeds op één krachtlijn. Een magnetische krachtlijn 
krijgt door de materie waar hij aldus aan gebonden is een identiteit 
en men kan gaan spreken van bewegende krachtlijnen en van verplaat- 
sing van fluxbuizen, wat in vacuüm of in een isolator zinloos is, 
maar ook bij de beschrijving van processen in supra-geleiders een 


bruikbaar concept is. 


3.2 Dissipatie 

Als men het eindige geleidingsvermogen in aanmerking 
neemt, ziet men dat de flux niet onbeperkt ingevangen kan zijn. 
Ingevangen flux betekane dat er stromen lopen; de EMK die deze 
stromen tegen de weerstand in drijft moet worden geleverd door een 
vermindering van de flux. Deze vervalt dientengevolge, en wel met 
een karakteristieke tijd tT, die we kunnen schatten aan de hand van 
het volgende voorbeeld. 

Laat een plasma in de X-richting begrensd zijn tussen 
x = 0 en x= TL en zich in de Y- en Z-richting oneindig uitstrekken, 


laat v = 0 en laat er een magnetisch veld in de Y-richting zijn: 
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Volgens VxBe= u is 
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NT. 


Hier is in de Maxwell-vergelijking de term oD/3t verwaarloosd tegen 


Ra | k > ne ei | 
J, wat in verband met nj = E + v Xx B neerkomt oo de veronderstelling 


t SSN [ vgl. (8) ] 


> + 
Volgens E = nj is 


_ _O X 
E‚ = En cos (7) s 
O 
> 
> > oB 
Volgens V x E = = == iS 

dt 
B_n JB 

Oo Ke Or Oo 
u_L? sin(5) = SE sin(7) 


zodat B, exponentieel afvalt met karakteristieke tijd 


u_LS 
Ee 
taiss — Tm 5 (5) 


Dit resultaat is hier afgeleid voor het geometrisch eenvoudigste 
geval om het proces ‘fysisch doorzichtig te maken. Mathematisch 


gezien loopt de afleiding langs de volgende lijnen: 
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Op zichzelf is dit resultaat niet nieuw; (5a) ligt ook 
ten grondslag aan een beschrijving van het bekende skin=effect, 
ten gevolge waarvan hoogfrequente stromen in een metaal slechts 


indringen tot een diepte die gegeven wordt door 
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Het begrip ingevangen of meer speciaal buitengesloten flux 
krijgt echter een bijzondere betekenis als we het in verband 
brengen met het begrip magnetische druk, waarover de volgende 
alinea's handelen. De tijd, tT, gegeven door (5) geeft dan een 
bovengrens voor de tijd gedurende welke een plasma in een der- 
gelijke eenvoudige configuratie door een magnetisch veld kan 
worden opgesloten en een ondergrens voor de frequentie waar- 


boven men van magnetohydrodynamisch gedrag kan spreken: 


u LZ . (6) 


3.3 Magnetische druk 
Een stationaire oplossing van (l) krijgen we door 


ET 0 te stellen. Dit geeft 


ÌxB-\p=0 … (7) 


Hieruit wordt j geëlimineerd door de Maxwell-vergelijking 


2, (s + 3 
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waarin de laatste Òerm wordt verwaarloosd, hetzij omdat we een 
stationaire toestand beschouwen, hetzij op grond van de in de 


magnetohydrodynamica gebruikelijke (en in 3.2 reeds gemaakte) 


onderstelling 1 


WS (8) 


Eo0 


Dan gaat (7) over in 


ar) dxE]xB-tp=0 





De laatste term is o.a. van belang bij kromming van de 
krachtlijnen, men kan hem dan beschouwen als een magnetische trek- 
spanning; in het bijzondere geval van rechte krachtlijnen waar- 
langs de veldsterkte niet verandert (zodat de gradiënt van de 


veldsterkte loodrecht op de krachtlijnen staat) wordt (9): 





B? _ 
v | Pp + 2e | = 0 … | (10) 


We zien hieruit dat het magnetische veld een druk ter grootte 
van B°/2u, op het plasma kan uitoefenen. Een voorbeeld waarop 
(10) van toepassing is, is de theta-pinch. Dit is een cilin- 
drische plasmakolom, opgesloten met behulp van een extern mag- 
netisch veld waarvan de krachtlijnen in axiale richting lopen. 
In het gebied waar de druk een (radiële) gradiënt heeft, loopt 
een azimutale stroom (vanwaar de naam Ô-pinch). Er bestaat ook 
een z-pinch, waarin de stroom axiaal loopt en het magnetische 
veld azimutaal. In de drukbalans hiervan kan men de laatste term 


van (9) in het algemeen niet verwaarlozen. 


3.4 Alfvên-golven 


Het feit dat het plasma en de magnetische krachtlijnen 
aan elkaar gekoppeld zijn geeft aanleiding tot golven, waarin 
de traagheid van het plasma en de trekspanning van het magne- 
tische veld een rol spelen. Beschouw een plasma met een sterk 
magnetisch veld B in de z-richting. Laat nu het plasma lood- 
recht op de krachtlijnen (b.v. in de x-richting) bewegen met een 


snelheid die een harmonische functie van z en tis (zie fig.2) 


en = v cos (kz —= wt) N (11) 


De hiervoor nodige kracht wordt geleverd door een stroom in de 


y=richting: 
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Fig. 2. Alfvên-golf. (In de tekst wordt B 


1 << B, ondersteld) 


of 
owv sin(kz —- wt) = Ì, X B ri 
zodat 


OP _ 
Je Bo wv sin(kz wt) (12) 


Met deze stroom hangt een verstoring van het magnetische veld, 


Wy 


‚’ samen. Hiervoor geldt 


md ae 
V Xx B SH 


of, daar we hebben aangenomen dat er geen x-afhankelijkheid 


ls: 


2 v cos (kz — wt) n (13) 
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Deze periodieke magnetische verstoring induceert een elektrisch 


veld volgens 


3B 
> > 
DS SE 
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Z dt 


waaruit, met behulp van (13), volgt: 


BP 2 
E = dj cos (kz = wt) … (14) 


Y B k? 


Tenslotte hebben we voor een ideaal geleidend plasma 


of 


hetgeen, met (11) en (14), oplevert dat 


= z V, . (15) 
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Deze grootheid, w/k, is de fasesnelheid van de golf. Men kan 

dit verifiëren door in het argument van de cosinus in (11) t/m 
(14) te substitderen zZz = wt/k. Deze keus van de snelheid van z 
maakt immers de cosinus, dus de fase van de golf, constant. 

Men noemt de grootheid B//u oP de Alfvén-snelheid naar de grond- 
legger van de magnetohydrodynamica, die in 1970 de Nobel-prijs 
ontving. 


Opgaaf: Laat zien dat de amplitude van de materieverplaatsing 


gelijk is aan de doorbuiging, &£, van de magnetische 


krachtlijnen (zie figuur 2). 
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Een wat meer formele en compacte behandeling van de 
Alfvén-golf verloopt als volgt: Men kiest de golfvector, kB en 
de snelheid, v si B Periodiciteit met exp i(k-r -wt) levert voor 
de differentiaal-operatoren: d/òt >» =iw en V > k. Verder geldt 
d/at = D/at +v-V = d/dt daar v 8 k. Tenslotte wordt “gelineariseerd"; 
alleen termen die van de Of of 1° orde zijn in de golf-amplitude 


worden meegenomen. 
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Bij de bovenstaande afleiding moeten we nog enkele 
kanttekeningen maken. Ten eerste mag men uit (15) niet afleiden 
dat de fasesnelheid naar oneindig gaat als p naar nul gaat. 


Als we de verplaatsingsstroom niet verwaarloosd hadden in (12) 


„zou de berekening niet (15) maar (15a) hebben opgeleverd (oefening!) ; 


de limiet van de fasesnelheid is daarom de lichtsnelheid. 

Ten tweede zij opgemerkt dat we een speciale en eenvoudige aanname 
hebben gemaakt voor het type golf (de mode); we hebben de verplaat- 
sing loodrecht op, en de voortplantingsrichting evenwijdig aan het 
constante veld B gekozen. De golf was dus transversaal en er trad 
geen compressie op. Maakt men de golf longitudinaal met voortplan- 
tingsrichting | B, (magneto-akoestische golf) dan verschijnt even- 
eens de geluidssnelheid in (15). In het meest algemene geval wor- 
den de formules vrij ingewikkeld. Ten derde hebben we de geleidbaar=- 
heid oneindig gekozen. Door dissipatie zal de golf — geheel analoog 
aan het voorbeeld dat tot (5) leidt - uitsterven met een kararteristie=- 
ke tijd T = HOL“ /n, waarin L = Kk L{oefening!). Eisen we dat T 

een flink aantal perioden bedraagt, dan moet wt >> 1 of 

WL u /N >> l, dus volgens (15): 


De dimensieloze grootheid R, Wordt het magnetisch Reynolds- 
getal genoemd. Men kan ook schrijven Re = V‚t/L. 

Men kan uit het bovenstaande concluderen dat het 
waarnemen van magnetohydrodynamische verschijnselen in het 
laboratorium niet gemakkelijk is, maar dat deze in de astro- 


fysica een dominerende rol spelen. Bijvoorbeeld: 


n = 1076 Ohm-m | n 
| (kwik) R, © 10 BL 
= 1,36 Xx 10* k À 
nj ge (B in Tesla, L in meter). 
Hos 4m x 10 ’ Henry/m 


Men moet dus ofwel een zeer sterk veld, ofwel een grote hoeveelheid 
kwik ter beschikking hebben. In gasontladingen kan de dichtheid 


bd 


b.v. een factor 10° ä 10° lager zijn, maar is de soortelijke 
weerstand heel vaak een factor 10% ä 10° hoger. Men moet daarom 
toch nog een flinke opstelling met een sterk geïoniseerd gas hebben 
om Alfvéên-golven te kunnen aantonen. Op astronomische schaal wordt, 


daar Te L?, R, in het algemeen zeer groot. 


eld 


4, Twee-vloeistoffen theorie 


In het vorige hoofdstuk werd het plasma beschouwd als een 





vloeistof (in de wijde betekenis die ook gassen omsluit) met elek- 

trisch geleidingsvermogen. Er zijn echter verschijnselen waarvoor 

dit theoretische model ontoereikend is, maar die wel goed beschreven 

| kunnen worden door het plasma te beschouwen als een mengsel van twee 

vloeistoffen, de ene met positieve en de andere met negatieve ladings- 

dichtheid, die ten opzichte van elkaar kunnen bewegen. Daarbij heeft 

; het ionen=-gas een hogere massadichtheid dan het elektronen-gas (voor 

een atomair waterstof-plasma al 1837 maal hoger), zodat met name 

| | voor hoogfrequente oscillaties in het plasma de ionen vaak als stil- 

staand kunnen worden beschouwd. In dit hoofdstuk wordt van deze be= 
nadering uitgegaan; het spreekt vanzelf dat men op de bereikte resul= 
taten correcties moet aanbrengen als men de snelheid van het ionen- 
gas wil verdisconteren. 

| | | De ladingsdichtheid van het elektronen=gas is "nae en de 

ME massadichtheid is nom’ als er in de bewegingsvergelijking alleen 

elektrische en magnetische krachten voorkomen, is de bewegingsver” 

| gelijking voor een volume-element van het gas identiek met die voor 

8 een vrij elektron. Hierop berusten de eerstvolgende drie paragrafen. 

| | Pas in Par. 4.3 zullen we verschijnselen ontmoeten waarbij ook de 

| druk van het elektronen=gas een rol speelt; dat die in 4.1 t/m 4.3 
niet in het spel is, komt omdat we daarin transversale golven be- 
schouwen die geen aanleiding geven tot V.v, dus ook niet tot druk- 

| gradiënten. Een andere reden om drukgradiënten te verwaarlozen zou 

| | kunnen zijn, dat de temperatuur laag is, wat volgens Pa = n.kTe de 

druk — en daarmee tevens de drukgradiënt — laag maakt; men spreekt 


| < dan van een koud-plasma theorie. 


4.1 Elektromagnet ische golven in plasma 


In deze paragraaf bespreken we de voortplanting van een 


\ 


gewone elektromagnetische golf door een ijl plasma. De afleiding 
is bv. van toepassing op radiogolven in de ionosfeer. Zoals reeds 
EO opgemerkt worden ionen-oscillaties verwaarloosd; de elektronen 


oscilleren volgens 
| Mm e= =-eË .… (16) 


13 


Deze beweging vertegenwoordigt een stroom 


+ > 
jj = = n_ev 8 
e” e 


In de Maxwell-vergelijking 
t DE 
je 
Vxbe Ho Ë "€ zE) 


krijgen we dus, behalve de verplaatsingsstroom die we in 
vacuüm hebben, ook nog een term j. We kunnen die twee samen- 
voegen door de bewegingsvergelijking (16) naar de tijd te dif- 


ferentiëren en (voor een harmonische golfbeweging met cirkel- 














2 
frequentie w) de operator „ te vervangen door wi. 
dt 
Dan krijgen we: 
Sed 
e mw? °t 
e 
zodat 
2 
en en 
mu? êt 
e 
waarmee 
2 
Ù ox  = Te DE 
7 kol So dt °° 
4 mw 


We zien hieruit dat voor kleine waarden van Ne de 
golfvoortplanting door het geïoniseerde gas overeenkomt met 
„ DN: : 
die in vacuum (waarin 5j=0) mits we Ee, Vervangen door 
B 5 noe“ /m w?, m.a.w. aan het plasma een relatieve diëlek- 
trische constante 


n e? 
e 
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toekennen. We kunnen de notatie vereenvoudigen door de plasma- 
frequentie in te voeren: 


ne 
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E Mm 
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(17) 
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ì (18) 


Zoals bekend krijgen we, door gebruik te maken van de 
andere Maxwell-vergelijking 


> ie 
Vv XxX Ez = 


een golfvergelijking waaruit voor vacuüm de lichtsnelheid 

ä = (oto) TS resulteert. Voor een diëlektrisch medium komt 
daar een factor er bij, m.a.w. de lichtsnelheid in ons plasma 
is 


Van Fm 


EM 7 8 
V l-=w 
p/® 


De brekingsindex, c/V is dus kleiner dan 1: 


EM’ 
Ce w‚ 
n=g=Ve= Vi E. 
E: w2 


4.2 Afscherming van de golf door het plasma 


Als de plasma-frequentie groter is dan de frequentie 
van de golf gaat in de Maxwell-vergelijkingen de term 5 over- 
heersen over 9D/at, waardoor e_ van teken omkeert, en de brekings- 
index imaginair wordt. De golfvergelijking 
92E 


> DE Ë B dE 2 À mn 2% 
Vv Xx V XE HoEofr ee of k°E Hoeoer® E 


krijgt dan, in plaats van de bekende periodieke oplossing, een 


exponentiële oplossing met 





k = WVE Ee = 5 


or Cc 
dus met een karakteristieke lengte 


de (mk) Vm san n 


wV|e,| _ wÄ-w? 
Ú p 
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De golf kan dan niet doordringen in het plasma, maar wordt af- 
geschermd in een grenslaag met een dikte d. Voor een zeer dicht 
plasma heeft d een asymptotische waarde 


d > : (19) 


Eg 
w 
P 


We hebben nu twee verschillende resultaten voor de 
afscherming van een E-M golf door een plasma. De berekening van 
het skin-effect voor een geleider met soortelijke weerstand n 


zegt — zoals we reeds zagen =— 


_ \ /2n 
a= V/a (20) 


zodat de vraag rijst: vanwaar het verschil? 

Het antwoord is opgesloten in vgl. (16), waarin we 
geen wrijvingsterm hebben opgenomen. Als echter de elektronen 
— door botsingen met ionen of neutralen —- een wrijvingskracht 
ondervinden, gaat de bewegingsvergelijking er als volgt uit- 


ziens: 





waarin ‚pn de tijd voorstelt, waarin bij afwezigheid van Ë de 


impuls door wrijving gedissipeerd wordt (To heet de botsings=- 


tijd, hi en de botsingsfrequentie). Is de botsingstijd zeer 


klein (wrs 1), dan domineert de wrijvingsterm rechts over de 


traagheidsterm links en krijgen we 





m_v 
a e 
ek + 2 = 0 
c 
of 
2 
5 ne En 
Js == neey = E « 
e me, 


De elektronensnelheid is dan in fase met E,‚, en wordt beschre=- 


ven met een soortelijke weerstand 
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e C _ Cc (21) 








e op 
We krijgen dan de resistieve skindiepte volgens (20). 
Zonder wrijving, maar met de traagheidsterm, krijgen 
we een stroom die 90° uit fase is met E en een skindiepte volgens 


(19), geldig als wt >> 1. In metalen is de botsingsfregquentie 


van de elektronen t.g.v. de hoge dichtheid van de ionen zeer 
hoog, zodat (20) meestal van toepassing is. Voor supragelei- 


ders geldt echter een indringing volgens (19). 


4.3 Faraday-effect 


In aanwezigheid van een magnetisch veld wordt de 
golfvoortplanting bijzonder ingewikkeld, daar de Lorentz-kracht 
op de oscillerende elektronen in aanmerking moet worden genomen. 
We nemen weer een eenvoudig geval, n.l. voortplanting van een 
transversale golf in de richting van het magnetische veld. Zo- 
als bekend kunnen we een vlakke golf ontbinden in twee circu- 
lair gepolariseerde componenten. Beschouw hiervan de rechts 


gepolariseerde component (w > 0): 


B, = E cos(wt-kz) ; E, = E sin(wt=kz) 


waarin X,Y,Z een rechtsdraaiend stelsel is waarvan de Z-as in 


| de richting van het constante veld B wijst. Kijkend in deze 


richting zien we de golf rechtsom draaien. 

We gaan nu weer gebruik maken van een notatie met complexe 
getallen. De conventie is nu echter anders dan in de voorafgaande 
paragrafen; we gaan een vector in het XY-vlak voorstellen als een 
complex getal, L=X+ iy en we voerende frequentie w in als | 
vector w = we. Dan geldt voor een in het XY-vlak met hoeksnelheid 


w roterende plaatsvector: 


Ke exp ilwt-kztd) 


> 
E 
©, 


> > > > 
Yr == liwr == WX r 
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waarbij w > 0 een rotatie in de zin van +X naar +Y voorstelt, 
dus rechtsom als men kijkt in de +Z-richting en w < 0 een draai- 
ing linksom. Verder mag k wel een complex getal zijn, maar dan 
blijft k een voortplantingsvector in de Z-richting en wordt geen 


vector in het XY-vlak. 


Met v_ = o, waarmee d/dt = d/dt, wordt de bewegingsvergelijking: 
Feöxie-S BB) -- SE d xt 
m o m ce 
e e 
waarin 
> _ EZ ee 
We En Bo (algemeen: Wo we B) (22) 


= 


de natuurlijke gyratiefregquentie van de elektronen om de magnetische 
krachtlijnen voorstelt. (oefening: ga na waarom in de bewegings- 
vergelijking de term vx Bolf verwaarloosd wordt) 


De bewegingsvergelijking kan worden herschreven als 


> > > e > 
(Ww -W _)Xr==-==E 
ce m 
e 
waaruit 
ee e > > 
F5 Tm (ww) aans: 
e ce 
en | 
n s noe’ ” Ee 
AT ae m_ (w-w_) à SC, 7 5 WW B Sg * 
e ce ce 


Nu levert Maxwell: 


- 
ri SO Kg dE 
OET vol ie: be) 


WW oe Z Z 
2 
Er 2E 
OO wlw-Woe) dt 
we 
zodat Ee = 2 ee | € 
w (w Woe) O 


A ik KWEENIE en nt rn erf AT edn enn in 


B Nik nn ee ks EE a il 


geldig voor rechts- en linksdraaiende golven als we positieve 
en negatieve waarden van w toelaten. Willen we w alleen als een 
positief getal zien, dan moeten we schrijven: 


à 





W 
En EE 
“rechts Ë www | Eo (23) 
ce 
en wi 
Bers It Mn = | E (24) 
links wW(wt ee Oo 


De twee circulair gepolariseerde componenten hebben blijkbaar 
een verschillende voortplantingssnelheid. Er treedt dan een 
draaiing van het polarisatievlak op (optische activiteit, 


Faraday-effect). Merk op dat de rechts-gepolariseerde golf een 


oneindige e (voortplantingssnelheid > 0) heeft voor w > Woe * 

Voorts dat € =0 voor wlw + w__) = w?, dus afhankelijk van 
gb 3 == CE P 

plasma-dichtheid (w_) en magnetisch veld (Woe): 


In een consequent doorgevoerde twee-vloeistoffen- 
theorie wordt ook de beweging van de ionen meegenomen. Er 
treden dan ook termen op met de resonantie-fregquenties 
Bi en w‚; van de ionen en de algemene formules voor golf- 
voortplanting in willekeurige richting worden bijzonder in- 


gewikkeld. 


4.4 Debye-lengte en plasmatrilling 


We besluiten deze beschouwingen over het twee- 
vloeistoffen-model met het introduceren van een karakteris- 


tieke lengte, de Debye -lengte, die in vele toepassingen een 


rol speelt. Laat de ionen en de elektronen in het plasma een 
ongestoorde, homogene dichtheid No hebben en laat de ionen=- 


dichtheid daarop gesuperponeerd een inhomogeniteit hebben 
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De elektronen zullen deze slechts zover opvullen, dat er nog 
een ladingsoverschot en dus een elektrisch veld overblijft, 
dat groot genoeg is om evenwicht te maken met de elektronen=- 
druk. 


Het ladingsoverschot is 


7 = e(n‚-n) = e(n, i°P,,e) 


waaruit een elektrisch veld resulteert: 


V.E = el 
E . a PE Ne 


en de drukbalans voor het (isotherm veronderstelde) elek- 
tronengas is: 


E = Vp = kr în = kT 3} 
ne = VPe = N= Hie . 
Samen: 
ne 
Vn = (n =- Nn. _) 
l,e Ee oKTe lt l,e 


zodat, als de verstoring een karakteristieke lengte L heeft: 





n — + = Li 
2 ik 2 
l,elL AB AD 
of 
L? 
Dive CS LZ + r dee 
waarin 
Ee KT 
M= 
D ne 7 


We vinden dus een karakteristieke lengte (de Debye-lengte) 
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(25) 


die door Debye werd ingevoerd in verband met elektrolyten. 


Als L >> Ap wordt Nie * Aji' de elektronen compen- 


seren dan de ruimtelading van de ionen. Ook als men externe 
velden aanlegt blijkt dat deze slechts een afstand in de orde 
van de Debye -lengte in het plasma kunnen doordringen. Een 
plasma waarvan de afmetingen groot zijn t.o.v. de Debye- 
lengte zal daarom steeds, behalve in de grenslaag, elektrisch 
neutraal zijn. 

De tijd die het plasma nodig heeft om deze neutra- 
liteit te bewerken volgt uit de bewegingsvergelijking: 


pn 
dv _ > 


Gecombineerd met 


> —_ OT _ B NN 
div(n ev) ket Ee) St (div E) 


levert deze (oefening!) een oscillatie met frequentie 





Ss (vgl. 17). 


Dit laat zien dat O5 ook de betekenis heeft van de natuur- 
lijke trilling (de plasmatrilling) die ontstaat op een 
plotselinge verstoring van het ladingsevenwicht. Als deze, 


b.v. door botsingen, is uitgestorven zal het evenwicht zijn 


hersteld. 
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5. Banen van individuele deeltjes 


In dit hoofdstuk zullen we nagaan hoe de individuele 
deeltjes reageren op elektrische en magnetische velden. Soms 
kan men hieruit zonder meer begrijpen hoe het plasma als ge- 
heel reageert, maar in andere gevallen stuit men eerst op, 
niet altijd even doorzichtige, statistisch-mechanische be-= 
rekeningen, In het eenvoudigste geval is het effect dat men 
beschouwt onafhankelijk van de thermische snelheid; men mag 
deze dan verwaarlozen. Effecten die evenredig zijn met de 
thermische snelheid, resp. energie, leveren na middeling over 
de snelheden een zelfde resultaat als men zou krijgen door te 
stellen dat alle deeltjes een middelbare snelheid, resp. 
energie, hebben. In dit geval, zowel als in het vorige, leiden 
beschouwingen gebaseerd op het gedrag van individuele deeltjes 
rechtstreeks tot de resultaten van de fluïdum-theorieën. 
Beschouwt men echter effecten die op een meer gecompliceerde 
manier van de thermische snelheid afhangen, zoals b.v. bij een 
resonantie van de thermische snelheid van een deeltje met de 
fasesnelheid van een golf (die aanleiding geeft tot een breuk 


met noemer v ze Vgore) dan moet men wel zijn toevlucht 


deeltje 


nemen tot de statistische theorie. 


5.l Constant en homogeen magnetisch veld 


De bewegingsvergelijking: 


Vv Vv 
X Y 
m Vv, = qB Me ë 
v 0 


Z 
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De Z-Component van de snelheid Vo ook vj genoemd, is dus 
constant, terwijl de X- en Y-componenten gekoppeld zijn. 
Eliminatie van één van de twee levert 


2 2 
… _ … (qB se 2e) 
He 2 vr resp. Vip En Le 


De oplossing hiervan is (het zijn er geen twee daar & wille- 
keurig isí) 


| Bt … Bt pi 
= VyCOS Ë dee + |: Vy = + Vasin |: de + 3 


waaruit blijkt, dat 


Eee vî (constant) 


en waarin de index | betekent: loodrecht op het magnetisch 
veld. De totale kinetische energie, km (vj + vs), is dus 
constant, wat niet verbaast, daar de Lorentz-kracht altijd 
loodrecht op de snelheid staat, dus geen arbeid verricht. 


We voerden reeds in de cyclotron-frequentie: 


en qB 
w 5 Oe (26) 
waarmee (twee notaties voor êên oplossing): 


v cos (wt + ò) cos (-w_t + &') 


= V‚ = Vi | 
: an s kan LJ 
‚2 sin(wt + 6) sin ( wat + ò') 


Dit is blijkbaar de natuurlijke gyratiefrequentie waarmee de 
deeltjes vrij langs de magnetische krachtlijnen kunnen spira=- 
liseren. Het teken in (26) geeft aan dat elektronen (met 


q = =e) rechtsom en positieve ionen linksom draaien als we in 


de richting van B kijken. 
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5.2 Diamagnetisme en drukbalans 


Zowel de gyratie van de ionen als die van de elek- 
tronen vertegenwoordigt volgens het bovenstaande een negatieve 
kringstroom. Het plasma is dus diamagnetisch. Elk deeltje 
heeft een magnetisch moment gelijk aan het produkt van de 
stroomsterkte en het omsloten oppervlak. De stroom is het aan- 


tal malen dat de lading per seconde de cirkel doorloopt: 


6 
Q 


Cc 
TT 


el 
Ii 


NO 
° 


De straal van de gyratiecirkel, For volgt uit 


Vv, = [w_ El 


waardoor het oppervlak van de baan 


mv? 


C 2 
Cc 


wordt,en we voor het magnetisch moment vinden: 


2 2 
AV MV 
U = == ge (27) 





Is de deeltjesdichtheid uniform en de fase willekeurig, dan vertegen= 
woordigt die gyratie geen nettostroomdichtheid. Waar een dichtheids- 
gradiënt aanwezig is, gaat echter wel een stroom lopen. Om deze 


te berekenen, nemen we even aan dat alle deeltjes dezelfde 


Y 


69 B, in het papier 


Vn 
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ie Eh he BBE En a ens ae aint ik bk nnn en. Sne hi md in lk dani hi nin 


gyratiestraal Fo hebben. Tot de stroom door het oppervlakte- 
element dxdz (zie Fig. 3) dragen dan deeltjes bij, waarvan de 


gyratiecentra in de gearceerde zone liggen: 


27 
wd 

, din 5 C 
JL | Dr n(x + r cos) dxdzd(r sing) 

$=0 

zodat 
2 wd 
ER an 
de = | DT [Go + r Cos 3e) r cos$dò : 
$=0 


Bij de integratie krijgen we alleen een van nul verschillende 
bijdrage van het produkt cos°$. Voor de berekening daarvan 


gebruiken we: 


27 2 
| cos? dè = | k(l + cos2b)db = TT 
0 0 


zodat 


ì se twart De an 
Jy ee dx '2B AX 





Gemiddeld over een thermische snelheidsverdeling is mví = 2KkT, 


zodat, wegens p = nkT, 
We zien dus, dat in dit geval 

ES X B — Vp = 0 
waarmee we aansluiting krijgen op (7). In een stationaire toe- 
stand wordt de. drukbalans blijkbaar mogelijk gemaakt door de 
diamagnetische kringstromen. OO 

In de bovenstaande afleiding hebben we dit laten zien 
voor het bijzondere geval dat de dichtheid een gradiënt heeft 


terwijl het magnetische veld en de temperatuur homogeen zijn. 


Een meer algemene afleiding van (7) zou vereisen dat we ook 
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de temperatuurgradiënt meenemen, waardoor de gemiddelde gyratie- 
“ straal een functie van de plaats wordt, en tevens de gradiënt in 
B, die twee gevolgen heeft,n.l. een gradiënt in de gyratiefre- 
quentie en een bijdrage tot de stroomdichtheid volgens vgl. (33). 
We laten deze afleiding achterwege; het bovenstaande voorbeeld 
van het effect van Vn illustreert voldoende hoe het verband is 
tussen macroscopische en microscopische beschouwingen in de 


plasmafysica. 


5,3 Constante en homogene elektrische en magnetische velden 


Heeft het elektrische veld een component in de rich- 


ting van B, dan worden de deeltjes in die richting onbeperkt 


versneld. Interessant is echter het geval É L B. Beschouw weer 


de bewegingsvergelijking 


> > > 
mv = q (E + v x B) … 


Een bijzondere oplossing (met v= 0) hiervan is: 
> > > 
E + v Xx B = 0 


of 


Vea (28) 


Voeren we nu in 





> > 
E Xx B 
B? 


. > > 
| Vv =v= 
| 


dan wordt de bewegingsvergelijking gereduceerd tot 


> > > 
m v* = qv x B 


ĳ h . f 
| m.a.w. het elektrische veld is weggetransformeerd en het pro- 
| bleem is gereduceerd tot dat van de vorige paragraaf. In ge 
kruiste É en B velden krijgen alle deeltjes - onafhankelijk 


van massa en lading - dus de "É x B driftsnelheid" volgens 


En > „ 
heden V‚y en vj volgens par. 5.1. Interpreteren we laatstge- 


| 

(28) en daarop gesuperponeerd de willekeurige constante snel- 
| 

| 
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noemde als thermische snelheden, dan zien we dat het plasma 
als geheel - koud of warm — loodrecht op E en B zijdelings 


wegloopt. Dit stemt overeen met (3). 
Merk op dat (28) alleen geldig is als de uitkomst vol- 


doet aan v, < c. Dit is geen zware eis aan B, want zelfs in een 
zeer zwak veld, als B = 107* T = l Gauss, mag E tot 30 kV mm! 
oplopen. Oefening: Ga na hoe een deeltje beweegt als niet aan 


bovenstaande eis voldaan is. 


5.4. Behoud van magnetisch moment 


Beschouw eerst een constant magnetisch veld, waarvan 


de krachtlijnen convergeren. Er is m.a.w. een V|IB| evenwijdig 


De ad 
aan B. Uit Fig. 4 is dan direct duidelijk dat de Lorentz-kracht, 


behalve een centripetale component, ook een axiale component 
heeft die in de richting van het zwakkere veld wijst en die we 
kunnen toeschrijven aan de radiële component van B (magnetische 


spiegel). Willen we het effect hiervan kwantitatief bezien, dan 


+) \ | 
Z > > ijk 
\ F= qVv Xx B 
B 
\ 
If © 
I I 
/ J 
V B 
„BI 


Fig. 4. Magnetische spiegel. 


moeten we de radiële component van B uitdrukken in VIB|, hetgeen 
kan met behulp van de wet V.B = 0, Volgens deze divergentie=- 
stelling moet de radiële flux door de cilinderwand van een plat 
cilindrisch doosje met straal r en hoogte Az, waarvan de sym- 
metrie-as evenwijdig is aan het magnetische veld, gecompenseerd 
worden door het verschil in de axiale fluxen door de platte 


eindwanden: 


Zi 





òB, 
2 Kms ame 
2nrAz<B > + TL 57 Äz 0 


of 


4 
SB, = 7 3 T5z 


waarin < > middeling over de omtrek aangeeft. Men moet eigenlijk ook 


ÒB, middelen, en wel over de eindvlakken van het cilindrische 


dz 
doosje, maar dit levert in eerste benadering niets anders dan 


B | | 
dat men 57 op de as moet nemen. Ter plaatse van de cyclotron= 


baan is, daar r_ = Vjre 8 


Cc C 
ien en ok ie 
r 2 AZ 2aB, da. 


De component van de Lorentz-kracht die met deze radiëele component 
van B samenhangt, zal een koppeling tussen longitudinale en 


transversale snelheid teweeg brengen: 


2 
mví òB 
Ee WSA 2B, 32 ° 

dv 


Deze vergelijking levert, na toepassen van jn = m TE en ver=- 
menigvuldigen met Vr gemiddeld over een cyclotronperiode: 


mv 3B 
_d Ee. UE: 
dt (5 mv) n 2B, zZz Vz * 


dw, We, 9B, , (29) 
dt B, dt | | 
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Nu wordt Len verwaarloosd in de veronderstelling dat 
de verandering van de gyratiestraal langzaam is ten opzichte 
van de gyratie zelf en t.o.v. de axiale snelheid. Dan moeten 
we -— daar de Lorentz-kracht geen arbeid verricht = op grond 


van behoud van energie verwachten dat 


W + W,‚ = constant 


z 0) 
zodat 
aw W. dB | | | 
wenn (30) 
Z 


Uit (29) zien we dat als het deeltje een magnetische spiegel 
nadert, zijn axiale snelheid afneemt, terwijl (30) laat zien 
dat daarbij de transversale snelheid wordt opgeslingerd. Uit 


(30) volgt tevens: 


dw dB 
et en Ee 
B, dt Wy dt 0 
of 
W 
d_ _Ò = 0 n 
de lef) - 


De grootheid tussen haakjes herkennen we uit (27), waaruit 


we zien dat 


al En constant , (31) 


Rl ee 


hetgeen geldig is bij Langzame verandering van de gyratie- 
straal. Merk op dat het deeltje (als We, # O0 is) slechts een 


beperkte diepte in de spiegel kan doordringen, immers 


W 
_ totaal 


W‚ < W ee 
Hm 


Ò totaal” 
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Behoud van magnetisch moment treedt ook op, als het 
veld niet in de plaats, maar in de tijd langzaam verandert. 


Dan wordt langs de gyratiecirkel een EMK geïnduceerd volgens 


2nr Er Sm Me DE 
SE 9B v, JB 
anton eat Baka. 5 
Ò 2 at 2w It  q2B At 


a 
ee 


| …_ dB _ 
waaruit, analoog aan het voorgaande (maar nu is rn 5E) 


QJ 


na middeling over een cyclotronperiode: 


W 
zodat wederom Ee constant is. 


Het inductieveld verricht wel arbeid, dus W‚, ver- 


andert terwijl W constant blijft. In plaats van de Ns Ò 
en z gebruikt men vaak, zoals we in het voorgaande reeds deden, 
indices L en //, die refereren aan de richting t.o.v. het mag- 
netische veld. De behoudswet voor het magnetische moment ligt 
ten grondslag aan het vasthouden van geladen deeltjes, resp. 
plasma in magnetische velden, b.v. bij thermonucleair onder=- 


zoek, in deeltjesversnellers en in de van Allen-gordels. 


> > 
5.5 \ Constant veld met VB L B 


Als de gyratiestraal een functie van de plaats is, 
gaan de helices die de deeltjes in een homogeen veld beschrij- 
van over in banen die - geprojecteerd op een vlak L B - er 
uitzien als de cycloïden van figuur 5. Kwalitatief kan men 
daaruit al zien dat er een zijdelingse snelheid loodrecht op B 
en VB zal ontstaan, waarvan het teken afhangt van de omloopzin, 


dus van het teken van de lading. 
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Fig. 5. Grad-B-drift, 


Ietwat schematisch aangeduid, verloopt de berekening 
van dit effect als volgt. We zoeken een oplossing van de bewe- 
gingsvergelijking van de vorm v= Vo + oe waarin Vv, de zijde- 


lingse drift moet worden en Vv. de daarop gesuperponeerde 


1 


cirkelbeweging. Met constante Vo wordt de bewegingsvergelij- 


king: 


8) (32) 


In het rechterlid zijn Vv en B constant en de andere twee, 


Ed > 


vj en EF» 


periodiek. Gaan we nu de periodieke termen sorteren, dan vin- 


periodiek met frequentie w. Het linkerlid is zuiver 


den we 


mt 


| 


@ 


nn 
m ad = Vv. Xx B + Vv x (Er ij B 
t 7 4 {Yi RR ed i 
Hieronder (33) blijkt, dat Ue evenredig is met v‚r;VB/B, 
zodat de twee termen in het rechterlid zich verhouden als 
(r‚VB/B)*. Door r‚ of VB klein te kiezen kunnen we de laat- 
ste term daarom willekeurig klein maken t.o.v. de andere, 


zodat overblijft 


d.i. de bekende vgl. voor een homogeen veld. De oplossing 
daarvan kennen we al; hij is periodiek met frequentie 


wW en Vv. = W L.. 
c 1 & 1 
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Om (32) kloppend te krijgen moeten we ook de niet=- 
periodieke termen sorteren. Schrijft men (oefening!) het produkt 
0 X ,.Ĳ)B in componenten, dan komt daarin van Vo zowel als 
van r,‚ een cosinus of sinus. Gemiddeld zijn cos? en sin? 


gelijk aan 5, dus de niet-periodieke termen uit (32) leveren 





Vaes + zv, r, VB = 
waardoor 
vr 
EN | | 


Het is van belang, op te merken dat ionen en elektronen in dit 
geval in tegengestelde richting bewegen, zodat er een stroom- 
dichtheid ontstaat ten gevolge van VB, evenals we eerder had- 
den ten gevolge van Vp (paragraaf 5.2). 

In de afleiding die volgt op (27) is, zoals aan het 
eind van 5.2 is opgemerkt, aangenomen dat B en T homogeen zijn 
en daarmee ook Wo en F+ Zouden we die afleiding correct wil- 
len doen met inachtneming van VB, dan zouden we o.m. ook de 
VB-drift volgens (33) mee moeten nemen om het juiste resul- 


taat te krijgen (oefening!). 


We kunnen tenslotte’ nog een verband aantonen tussen 
de effecten die veroorzaakt worden door V‚|B| (deze paragraaf) 
en VIB (paragraaf 5.4). 

P In 8 5.4 hebben we gezien dat Vv [B| // B een kracht 
in de z=-richting oplevert, die gegeven wordt door 


of OO (34) 
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De driftsnelheid volgens (33) kan eenvoudig worden herleid 





tot 
v? mv; Mar 
yv = VB = B VB_ = aB VB, 
o 2B_W Z 2aB, Z aB, 


of, als men de richting van Vv in aanmerking neemt: 


E 


CN: (35) 
Vo zz Bx V [Bl . 


Q 


Nu leert een afleiding, geheel analoog aan die van (28), dat 

> El: 5 B | n 

een kracht F,‚ een driftsnelheid - oplevert. Blijkbaar 
qB? 

kan (35) ook worden verkregen door een kracht 


Va |B (36) 


in te voeren, analoog aan (34). 


Aanbevolen boeken: 


Hoofdstukken over plasmafysica in 
Reitz en Milford, Foundations of Electromagnetic Theory 
(korte inleiding). 
J.D. Jackson, Classical Electrodynamics (uitvoeriger). 
Verder: 
DJ. Rose and M. Clark Jr., Plasmas and Controlled Fusion. 
G. Schmidt, Physics of High Temperature Plasma. 
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Lijst van symbolen 


tijt Mt Wy 


Hd 


m 


=, 


HD > 
Es 


Fh 


A 


Q 


« 
G- 
« 


magnetische veldsterkte (inductie) 
elektrische veldsterkte 

kracht 

het imaginaire deel van 
karakteristieke lengte 
oppervlakte 

magnetisch Reynolds-getal 
temperatuur 

Alfvên-snelheid 


voortplantingssnelheid E-M golf 
fasesnelheid 


kinetische energie 

lichtsnelheid 

index voor "cyclotron" in Far Wo 
index voor "collision" in Vor Te 
indringdiepte 

elektronlkading, 1,6 X 1019 C 
index voor elektron, resp. ion 
/-1 

elektrische stroom 
stroomdichtheid 

constante van Boltzmann 


golfvector, golfgetal 


l 
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massa van een deeltje (m = 9,11 Xx 10 
dichtheid van deeltjes 
druk 
lading van deeltje 
cyclotronstraal 
tijd 
snelheid 
cilindercoördinaten 
cartesische coördinaten 
index voor evenwijdig, resp. loodrecht 
diëlektrische constante 
soortelijke weerstand 
Debyelengte 
magnetische permeabiliteit 
magnetisch moment 
botsingsfrequentie 
34 
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verplaatsing 

massa per volume-eenheid 
lading per volume-eenheid 
botsingstijd 

magnetische flux 
zwaartekracht-potentiaal 
cirkelfregquentie 
cyclotronfregquentie 


plasmafreguentie 
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